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■ Abstract 

' The production mechanism of highly excited nuclei in the Fermi energy domain is investi- 

, gated. A phcnomenological approach, based on the exciton model, is used for the description 

of pre-equilibrium emission. A model of deep inelastic transfer is employed for the peripheral 
collisions in the post-pre-equilibrium stage. An approach to describe more central collisions is 
' proposed. A geometric overlap formula is employed in a way suitable for given energy domain. 

\Q , A simple geometric approach describing the interaction of participant and spectator zones is 

used to determine the incomplete fusion channel. Excitation energies of both fragments are 
determined. Results of the calculation are compared to available experimental data and an 
overall satisfactorily agreement is obtained. The model's ability to describe the production of 
, the hot nuclei can be employed in the study of multifragmcntation and/or in the production of 

rare beams. 
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A detailed knowledge of the production mechanism of hot nuclei is desirable for the under- 
standing of multifragmentation. The peripheral collisions of heavy ions in the Fermi energy domain 
demonstrate that the nucleon exchange is a dominating production mechanism of the highly excited 
quasiprojectiles. In more violent collisions, the processes leading to the mid-velocity emission start 
to play an important role. The influence of the pre-equilibrium emission and fragmentation-like 
rS \ processes becomes an important issue. There exists a considerable amount of experimental data 

from the reactions of heavy ions in the Fermi energy domain ( for review see e.g. refs. |ffl, 0) ). 
Various transport and molecular dynamics models are generally used for a theoretical description 
( for review see e.g. ref. || and references therein ). In this article, a framework based on simple 
phenomenological assumptions will be presented and its ability to describe the experimental data 
will be demonstrated. 

The pre-equilibrium emission ( PE ) is a process where fast particles are emitted prior to 
the equilibration of the system. Emission of pre-equilibrium particles in reactions induced by 
nucleons and light particles was theoretically explained using the exciton model @ or using the 
master equations || . For reactions induced by heavy ion beams a model of nucleon exchange was 
developed ||. In the present work, we use a phenomenological description j7j based on similar 
assumptions as the exciton model. The probability of pre-equilibrium emission is evaluated using 
the formula 

(n/neq— I) 2 

P P re{n/n eq ) = 1 - e (1) 
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for n < n eq and is assumed zero for n > n eq , where n is the number of excitons at given stage and 
n eq is the the number of excitons in the equilibrium configuration for given excitation energy. The 
a is a free parameter. The basic assumption leading to eq. ( [l] ) is the dependence of P pre on the 
ratio n/n eq as suggested in ref. (8|. An initial exciton number is equal to the mass number of the 
beam. The equilibrium number of excitons is calculated according to the formula ||] 

n eq = 2gTln2 (2) 

where g is the one particle level density at the Fermi energy and T is the nuclear temperature 
determined as T 2 = U/d, where a is the level density parameter and U is the excitation energy. 
At every emission step, a random number between zero and one is generated. If the random 
number is smaller than P pre a pre-equilibrium particle is emitted. The emission of neutron, proton 
and a-particle is considered and the Weisskopf-Ewing emission widths are used. The Maxwellian 
spectrum of kinetic energy with the apparent temperature |J 

T app = l^(E P - Vc)] 1 ' 2 (3) 

is assumed, where Ap is the projectile mass number, Ep is the projectile energy and Vc is the 
Coulomb barrier. The emission angle is determined according to the formula [1C] 



where A9 = k ^ CN , Rcn is the radius of the compound nucleus and k is the wave number of the 
emitted particle. After emission, the exciton number is increased by a value obtained using the 
formula 

An = A pre - — (5) 

Prad 

where A pre is the mass of emitted particle, (3 ra( i is the radial velocity in the contact configuration 
at a given angular momentum and k is a free parameter. If no pre-equilibrium emission occurs at 
a given emission stage, the pre-equilibrium stage is finished. 

As a next step, the interaction of the projectile and target is considered. As shown in several 
works [11, 12 1, the model of deep inelastic transfer ( DIT ) describes well the peripheral collisions 



where the relative motion is mostly tangential. With decreasing angular momentum the radial 
motion becomes more intense and violent scenarios should be considered. At low energies several 
MeV/nucleon above the Coulomb barrier an incomplete fusion was observed by detecting evapora- 
tion residues |l0| , |l3| , 14] or fission fragments [jl5). The framework of angular momentum windows 



was introduced by Wilczynski et al. At very high projectile energies the abrasion-ablation 

model [O] is used widely to describe the fragmentation phenomena. In the proposed description, 
different scenarios are employed depending on angular momentum. At large angular momenta, the 
model of deep inelastic transfer is used. In the more central collisions, where stationary di-nuclear 
configuration cannot be created, the framework of a geometric overlap is used. In the near-central 
collisions a model of Wilczynski 0] is employed. 

Since pre-equilibrium emission occurs prior to the fragmentation stage, it is necessary to recon- 
struct the post-pre-equilibrium projectile-target configuration. It is assumed, according e.g. to the 
conclusions of work || that pre-equilibrium particles are mostly emitted from the ( usually light 
) projectile and propagate through the target. The emitted mass and charge are subtracted from 
the projectile. The excitation energy of the target is set equal to the sum of kinetic energies of 
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the emitted pre-equilibrium particles, in accordance to the assumption that one nucleon-nucleon 
collision occurs during the propagation of the particle through the target. 

For every event the Monte Carlo DIT code of Tassan-Got |11] is used. In the case when di- 



nuclear configuration is created, deep inelastic transfer takes place and excited quasi-projectile and 
quasi-target are created. In the cases, where overlap of nuclei is too deep, the more violent scenario 
( here arbitrarily denoted as realistic geometric fragmentation - RGF ) is considered. The geomtric 
overlap formula of the abrasion- ablation model [17| is used. It is clearly unrealistic to assume 
that the projectile propagates along straight line determined by the asymptotic value of impact 
parameter. Instead, a minimum distance between projectile and target in the Coulomb scattering 
of two point-like charges is used. Thus, one participant and one or two spectator zones are created. 



The charges of the spectators are determined according to the combinatorial formula 1|, 19] 



p(z iS ) = $ S )(£; S )/C&) (6) 

where i=P(T) for the projectile ( target ), Ap( T ), Zp( T ), N P m are the mass number, charge and 
neutron number of the projectile ( target ) and Ap(p) S , Z P ^ S , Np( T )s are mass, charge and the 
neutron number of the projectile ( target ) spectator. The charge of participant zone is determined 
assuming charge conservation. 

In the Fermi energy domain one can assume that the participant zone can be captured by either 
projectile or target spectator zone. In order to decide between these two, the volume occupied by 
the neighboring nucleons within the reach of nuclear interaction ( 1 fm ) is determined in both 
spectators. The volume is approximated by a 1 fm thick slice of the sphere. The number of 
neighboring nucleons ( A^s ) is then determined using a normal distribution centered at the value 
exactly corresponding to the volume with the standard deviation equal to \AAjvs- The numbers of 
neighboring nucleons are compared and the participant zone is captured by the spectator with more 
neighboring nucleons. Such a procedure reflects the saturation of the nuclear force. The capturing 
spectator and the captured participant zone form a very hot fragment. The other spectator is much 
colder. In this case, the excitation energy is determined using the formula 

E* s = xA NS (^--Vc>)^ (7) 

where Ep/ and Ap/ are the kinetic energy and the mass number of the effective projectile after pre- 
equilibrium emission, Vc is the Coulomb barrier, < s > is the mean effective path of the nucleon 
along the spectator trajectory within the slice, A is the mean free path between the nucleon-nucleon 
collisions in the nucleus ( a value 6 fm is chosen ) and x is a random number between zero and one. 

The energy and the angle of the spectator is determined using the formula of Matsuoka et al. 
[20] based on the Serber approximation. The formula reads 



d 2 a (E a E b )V* 



dE a dti a (2^lB p > + 2m 2 a E P i/m P i + 2m a E a - A(mlEp>E a /m p/ ) 1 / 2 cos9) 2 

where it is the fragment a which flies away and the fragment b which fuses with the other nucleus, 
E a and Et, are their kinetic energies, Bpi is the binding energy of a and b in P, [i is the reduced 
mass of the system a + b , mp>, m a , rrib are the masses of P',a,b and 9 is the emission angle of a with 
respect to the direction of P'. The kinetic energy and angle of a hot fragment are determined using 
energy and momentum conservation laws. An intrinsic angular momentum is calculated using a 
mean radial distance and momentum of the participant zone relative to the capturing spectator in 
the configuration of minimum distance between the projectile and target. 

In order to describe also the inverse kinematics, namely when the projectile is heavier than the 
target, the system is transformed into the inverse frame where the projectile becomes a target and 
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vice versa. Then the calculation proceeds as described above and the final kinematic properties of 
the reaction products are obtained after a re-transformation into the lab frame. 
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Figure 1: Experimental ( symbols ) and calculated ( lines ) mean multiplicities of pre-equilibrium neutrons as 
a function of kinetic energy loss of the projectile-like fragment. Solid squares - experimental multiplicities measured 
in reaction of 35 MeV/nucleon 48 Ca beam with 112 Sn target, open squares - ditto for 35 MeV/nucleon 40 Ca beam, 
solid line - calculated multiplicities in reaction of 35 MeV/nucleon 48 Ca beam with 112 Sn target, dashed line - ditto 
for 35 MeV/nucleon 40 Ca beam. 

The model of pre-equilibrium emission was compared to the results of work |2l]] where a multi- 
plicity of the pre-equilibrium particles was determined in coincidence with projectile-like fragments 
( PLFs ) in the reactions of Ca beams with 112 Sn target at 35 MeV/nucleon. In Fig. 1 are given the 
values of pre-equilibrium neutron multiplicity in reaction 40 ' 48 Ca-|- 112 Sn for several bins of kinetic 
energy loss of the projectile-like fragment. The solid ( open ) squares represent the results of work 
[21] and the lines represent the results of the calculation. The agreement is quite good. The pa- 
rameters cj=0.25 and k=0.3 were used in the calculation. The same values of a and k were tested 
in other reactions and lead to results which track well with the results of experimental works where 
multiplicities of pre-equilibrium particles were determined in coincidence with heavy residues or 
fission fragments |22|, ^] or in coincidence with the quasi-projectile [p^|. 



In the recent experimental work [ 23 1 a linear correlation between the primary mass of the 
projectile-like fragment and the net mass loss due to the de-excitation was reported in the nearly 
symmetric reactions of 93 Nb with 116 Sn at 25 MeV/nucleon in both normal and inverse kinematics 
for different dissipation bins. The net mass loss increases with the primary mass of the projectile- 
like fragment. With increasing dissipation this trend occurs in the still broader range of primary 
masses. Since the net mass loss is possibly correlated to the excitation energy of the hot primary 
projectile-like nucleus, one may expect similar trend also for excitation energy. Such a trend is 
a possible signal of the breakdown of the concept of deep inelastic transfer. In Fig. 2 is given a 
calculated correlation between the excitation energy and the mass of the hot projectile-like nucleus 
for different bins of kinetic energy. One can see that the calculation follows the experimental 
trend. At masses close to the beam the deep inelastic transfer takes place but the range of primary 
masses is quite narrow. To achieve larger mass change a more violent collision should occur. When 
the target strips a part of the projectile, the projectile-like fragment remains relatively cold. Hot 
projectile-like fragments are produced if a part of the target is picked-up by the projectile. 

The production of heavy residues was studied recently by Skulski et al. 0] in the reaction 
86 Kr+ 197 Au at projectile energy 35 MeV/nucleon. A spectrum of kinetic energies of the target-like 
fragments ( TLFs ) was measured at angles 9 - 46 ° in coincidence with projectile-like fragments 
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Figure 2: Calculated correlation between the primary mass and the excitation energy of a projectile-like fragment 
in the reaction 93 Nb with 116 Sn at 25 MeV/nuc 
PE+DIT events, grey squares - PE+RGF events. 



in the reaction 93 Nb with 116 Sn at 25 MeV/nucleon. Four different bins of E PLF are defined. Black squares 



with Z > 25 at angles 2 - 8 °. The spectrum exhibits two humps, one of them is the low energy 
component at energies below 40 MeV and the second one is located between 50 and 100 MeV. Fig. 
3a shows the calculated spectrum of kinetic energies for the same reaction using identical angle 
cuts and charge threshold. One can see that similar range of the kinetic energies is covered. The 
double- humped structure is not as clearly pronounced as in the experiment. As may be seen in Fig. 
3b, the double-humped shape can be explained by a correlation between the kinetic energy and the 
angle. The low energy particles are emitted at large angles and an additional energy losses or re- 
scattering can be caused by interaction with the target. The high energy part is possibly influenced 
by decreasing mass and charge and increasing excitation energy of the coincident projectile-like 
fragment what leads to the decrease of detection probability because of the charge threshold. 
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Figure 3: Calculated kinetic energy spectrum of target-like fragments in reaction 86 Kr+ 197 Au at projectile energy 
35 MeV/nucleon and the correlation of Otlf and E' rLF for the same reaction. For details see text. 
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The production of heavy residues in inverse kinematics was measured recently by Souliotis et al. 
pq| in the reaction 197 Au+ na *Ti at projectile energy 20 MeV/nucleon. Fig. 4 shows the measured 
yields of heavy residues at the forward angles as a function of A and Z. The solid line represents 
the calculated centroids of Z for a given residue mass. The code GEMINI |2(|] was used for the 
de-excitation stage. One can see that the calculation follows the experimental trend quite well. In 
this case the participant zone sticks almost exclusively to the heavy projectile and the incomplete 
fusion leads to the production of neutron-deficient residues with masses close to the beam. 

The measurement of the production of intermediate mass fragments ( IMFs ) in symmetric col- 
lisions 58 Fe, 58 Ni+ 58 Fe, 58 Ni at 30 MeV/u [§7| determined three different sources of IMFs, the mod- 
erately excited projectile(target)-like source at velocities near the projectile (target) and the highly 
excited source at velocities near the center of mass velocity. Fig. 5 shows a correlation between the 
excitation energy and the velocity of reaction products in the lab frame. Both projectile-like and 
target-like nuclei are included in the plot. One can identify three sources analogous to the ones seen 
in experiment. The projectile- and target-like sources are moderately excited. The third source 
with the average velocity close to the mid-velocity is highly excited. Both the PLFs and TLFs 
are included in this source. The difference of the isospin of IMFs from different sources which was 
experimentally observed may possibly be explained by different excitation energies of the sources. 
As follows from refs. [28, 29], the isospin dependences of fragment yield ratios are more flat at high 
excitation energies and the sensitivity to isospin of the source is weaker. 
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Figure 4: Measured yields |g5[ of heavy residues at the forward angles in the reaction 197 Au(20 

MeV/nucleon)+ nat Ti as a function of A and Z. Z is expressed relative to the line of /3-stability. Solid line - cal- 
culated centroids of the fragment charge for given residue mass ( the code GEMINI [^| was used for de-excitation ). 

The comparison to the broad range of experimental observables measured in various reactions 
in the Fermi energy domain appears to imply that the present approach describes correctly the 
processes leading to the production of excited projectile-like and target-like nuclei in the range 
between 20 and 50 MeV/nucleon. With increasing projectile energy, the production of three-body 
events and the compression phenomena start to play an important role. As a primary candidates 
for three-body events can be considered hot fragments where the relative motion of the participant 
zone leads to values of intrinsic angular momenta above the critical angular momentum for fusion. 
The compression phenomena can primarily take place in the events where the fragment with the 
mass close to the compound nucleus is produced with low intrinsic angular momentum. 

In conclusion, the approach presented in the present work appears to be a suitable tool to 
describe the production of hot nuclei which later undergo multifragmentation. Furthermore, pro- 
duction of rare beams in the Fermi energy domain can be addressed using this approach ( for 




6 




-0.05 0.05 0.1 0.15 0.2 0.25 0.3 



Figure 5: Calculated correlation of the excitation energy and velocity of the hot nuclei produced in the reaction 
58 Ni+ 58 Ni at 30 MeV/nucleon. 

instance, experimental studies are underway at the Cyclotron Institute of Texas A&M University 

H ). 
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